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Dynamics and correlation length scales of a glass-forming liquid in quiescent and 

sheared conditions 
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We numerically study dynamics and correlation length scales of a colloidal liquid in both quiescent 
and sheared conditions to further understand the origin of slow dynamics and dynamic heterogeneity 
in glass-forming systems. The simulation is performed in a weakly frustrated two-dimensional liquid, 
where locally preferred order is allowed to develop with increasing density. The four-point density 
correlations and bond-orientation correlations, which have been frequently used to capture dynamic 
and static length scales £ in a quiescent condition, can be readily extended to a system under 
steady shear in this case. In the absence of shear, we confirmed the previous findings that the 
dynamic slowing down accompanies the development of dynamic heterogeneity. The dynamic and 
static length scales increase with a-relaxation time t q as power-law £ ~ r£ with fj, > 0. In the 
presence of shear, both viscosity and r a have power-law dependence on shear rate in the marked 
shear thinning regime. However, dependence of correlation lengths cannot be described by power 
laws in the same regime. Furthermore, the relation ^ ~ t„ between length scales and dynamics 
holds for not too strong shear where thermal fluctuations and external forces are both important 
in determining the properties of dense liquids. Thus, our results demonstrate a link between slow 
dynamics and structure in glass-forming liquids even under nonequilibrium conditions. 

PACS numbers: 61.20.Ja, 61.20.Lc, 64.70.P- 



I. INTRODUCTION 



Glass can be formed in a variety of molecular and col- 
loidal systems, but the nature of the glass transition re- 
mains elusive despite intensive research [D-Ql- It is com- 
monly observed that time scales, such as structural re- 
laxation time, dramatically increase on the approach to 
the glass transition point with little change in the static 
structure of the liquid. This is the reason why the re- 
search on the glass transition for a long time has been 
focused on time scales much more than length scales. 
In the past decade, however, both experiments and sim- 
ulations have revealed that the relaxation process not 
only drastically slows down but also becomes progres- 
sively more heterogeneous as a liquid approaches its glass 
transition point [a-Hlf. The collective nature of the dy- 
namics of the supercooled liquids can indeed date back 

and has been 
II . The phe- 



to the seminal idea of Adam and Gibbs |12| 
well reviewed in recent publications [H, Il3r 
nomenon of growing dynamic heterogeneity approach- 
ing the glass transition point also suggests the existence 
of a growin g d ynamic length scale in glass-forming sys- 
tems 

BEL Ell- As 

a matter of fact, the issue of charac- 
teristic length scales has become one of the main themes 
of the study on the glass transition Q ■ 

To further understand the origin of the dynamic het- 
erogeneity and explore whether there exist structural sig- 
natures connected to the slow dynamics and dynamic 
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heterogeneity in viscous liquids, several ideas have been 
put forward. One of the ideas is based on the effect of 
frustration on crystallization in vitrification. To continu- 
ously control frustration, several ways are frequently used 
in the study fl8rl20j . One way is to focus on the polydis- 
perse colloidal liquid where frustration can be continu- 
ously changed by varying polydispersity [l9[ (as inspired 
by the study of crystal nuclcation on polydispcrse col- 
loids [H|) and another way is to study liquids in nega- 
tively curved space where frustration can be controlled 
by varying space curvature (2(| (as motivated from the 
frustration-based approach of the glass transition fl8j). 
In both cases, if the frustration is not very strong, some 
local order reminiscent of the crystal structure can de- 
velop with increasing density or decreasing temperature. 
Then, one can readily detect some static correlations 
(usually by bond-orientation correlation functions) ap- 
proaching the glass transition and investigate whether 
these structural signatures are connected to the dynam- 
ical properties. Although different interpretations were 
given for the change of correlation length scales approach- 
ing the glass transition point [l9|, [20||, the structural 
signature of slow dynamics and dynamic heterogeneity 
does exist at least in these model glass-forming liquids. 
In particular, a direct relation between dynamic hetero- 
geneity and the so-called medium-range crystalline order 
(MRCO, i.e., long-lived clusters of particles with high 
structural order) has been revealed in colloidal disper- 
sions [lj|. Another more general method to detect the 
nontrivial static spatial correlations has been proposed 
recently, which is motivated from the random first or- 
der transition (RFOT) theory [22|, [23| (see recent re- 
views d, HU for details on this theory) . Within this the- 
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ory, the so-called point-to-set correlations |25l - l29{ are be- 
lieved to capture nontrivial structures of viscous liquids. 
The idea is to first freeze the position of a set of parti- 
cles in an equilibrium configuration and then let the sys- 
tem evolve in the presence of the constraint and measure 
how the position of the remaining particles is affected. 
Thus, the point-to-set length is a measure of the spatial 
extent over which the effect of equilibrium amorphous 
boundary conditions propagates. Numerical simulations 
have confirmed qualitatively the growth of the point-to- 
set correlations |25| and its connection to the dynamic 
correlations has also been explored recently (2(| [27J ■ Be- 
sides, the structural features of slow dynamics in viscous 
liquids have also been investigated in Refs. [30l - f36| . 

In recent years, the study of supercooled liquids and 
glasses under nonequilibrium conditions (e.g., under 
shear) has also attracted considerable attention, with em- 
phasis on their rheological properties and the establish- 
ment of the concept of effective temperature j37l - l4l| . The 
rhcology of soft glassy materials not only provide new in- 
sight into the nature of the glass transition but also has 
its own value. One common observation is the marked 
shear-thinning behavior in glass-forming liquids. In par- 
ticular, the shear viscosity rj is found to decrease with in- 
creasing the shear rate 7 as 77 ~ with v < 1 . Another 
intriguing finding is that soft glassy materials commonly 
display the shear banding phenomenon [39j . i.e., a sys- 
tem under strong shear spontaneously "phase separates" 
between a flowing region supporting the shear and an im- 
mobile region with no flow. Recently, the question of how 
dynamic hetero geneity is influenced by the flow has also 
been addressed j42T - |48| . e.g., numerical simulations HU 
revealed the existence of large scale heterogeneities under 
shear flow and its connection to the non-Newtonian be- 
havior of a supercooled liquid has also been explored, and 
power-law dependence of the intensity of the dynamic 
heterogeneity on shear rate has been observed [44| - |48| . 
However, it is still not clear how dynamic heterogeneity, 
MRCO and the associated length scales respond to shear. 

In this work, we focus on slow dynamics and dynamic 
heterogeneity of a colloidal glass-forming liquid and their 
structural features in both quiescent and sheared condi- 
tions. We use a weakly frustrated two-dimensional (2D) 
liquid as our model system, where crystallization can be 
avoided and locally preferred order is allowed to develop 
with increasing density. The four-point density correla- 
tions and bond-orientation correlations, which have been 
frequently used to capture the dynamic and static length 
scales £ in a quiescent condition, can be readily extended 
to a system under steady shear in this case. We con- 
firmed the previous findings in the absence of shear that 
slow dynamics accompanies the development of dynamic 
heterogeneity. Both dynamic and static length scales in- 
crease with a-relaxation time r Q as power-law £ ~ 
with \i > 0. In the presence of shear, both rj and r Q have 
power-law dependence on 7 in the marked shear thinning 
regime. However, dependence of the dynamic and static 
correlation lengths on 7 cannot be described by power 



laws in the same regime. We further find that the relation 
£ ~ T a between length scales and dynamics holds for not 
too strong shear where thermal fluctuations and applied 
forces are both important in determining the properties 
of dense liquids. The paper is organized as follows: In 
Sec. II, we show the model and methods used in this 
study, and then we present results and discuss structural 
signatures of slows dynamics and dynamic heterogeneity 
in the absence of shear in Sec. III. Results under steady 
shear are given in Sec. IV. Finally, Sec. V summarizes 
our results. 



II. MODEL AND METHODS 

The model we focus here is a 2D liquid, in which 
the particles interact via the Weeks-Chandler- Andersen 
(WCA) potential H|: 

a ., (r) = Ni^k/r) 12 - {o, jk /rf + 1/4], for r < 2^a jk 
J (0, otherwise, 

(1) 

where e is the depth of the potential well, r is the dis- 
tance between two particles, and Ujk = (cj + <Xfc)/2 with 
<jj the diameter of particle j. As the WCA potential 
is short-ranged repulsive, a WCA liquid at low tempera- 
tures behaves as a colloidal liquid. Although recent simu- 
lations [5CJ, l5l| have shown that the longer-ranged attrac- 
tive interactions can affect the dynamics of viscous liquids 
in quantitatively and qualitatively ways, we expect that 
the conclusions in our work will not be changed since the 
purely repulsive model can also capture the essential fea- 
tures of a typical glass-forming liquid. To introduce weak 
frustration on crystallization in vitrification, the parti- 
cle diameters were chosen equidistantly from the range 
0.8 — 1.2 with an interval of 0.01, then the pol ydispcrsity 
for the system is A = y/{< <r 2 > — < a > 2 )/ < cr >= 
11.98%, where < x > is the average of variable Xj among 
all the particles. Such polydispcrsity can avoid crystal- 
lization in the system and allow static structural order to 
develop with increasing density. In principle, Brownian 
dynamics simulation is a more faithful method to study 
colloids. However, previous work has suggested that at 
long times, when particles have collided many times, 
Brownian and Newtonian systems show qualitatively the 
same phenomenology [52j . We employed molecular dy- 
namics simulation in the NVT ensemble, where Newton's 
equations of motion were integrated for N — 1000 parti- 
cles with the velocity form of the Verlet algorithm under 
periodic boundary conditions and the temperature T was 
maintained by the Nose- Hoover thermostat [HH . All the 
particles have the same mass m. Length and time were 
expressed in units of < a > and yfm < a > 2 Je. The 
time step is At = 0.002. For convenience, the tempera- 
ture T in our simulation is fixed at ksT/e = 0.025 with 
ks the Boltzmann's constant. We use the area fraction 
4> = (l/V) x T,f=in(vj/2) 2 (V = L 2 with L the box 
dimension) as a control parameter. In this work, we did 
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not make a correction to the softness of the WCA po- 
tential and map a WCA system to a hard-sphere system 
by using the effective diameter since we only want to fo- 
cus the glassy behavior as the density varies. In athermal 
colloids, plays the same role as temperature in molec- 
ular liquids. The area fraction <fi in this work covers the 
range from 0.62 to 0.65 (L = 35.01 ~ 35.85). In order to 
check the possible finite size effects, we also performed 
simulations for a system with N = 4000 in the absence 
of shear, and we found the qualitative agreement with 
the results from the smaller system. Thus, we only show 
the results with N = 1000 in this paper. At each 4>, the 
system was equilibrated for 5 x 10 7 time steps first before 
collecting data and 8 independent runs were performed. 

In order to study the system under steady shear, we 
used the nonequilibrium molecular dynamics (NEMD) 
SLLOD algorithm 0, The SLLOD equations of 

motion are given by: 



dt 



m 



Tj ■ Vv, 



dp, 

dt 



(2) 



where Pj and are the peculiar momentum and labora- 
tory position of particle j, and Vv is the gradient of the 
streaming velocity v. Combining with the Lees-Edwards 
boundary conditions (56j . these equations of motion can 
generate homogeneous flow, which means that no insta- 
bility such as shear banding (39j occurs in our simula- 
tions. In the case of planar Couette shear flow, v has 
only one non-zero component, which is given by 71/x, 
where 7 is the shear rate and x the unit vector along the 
x axis. Here the velocity gradient is in the y direction, 
and the streaming velocity in the x direction. The shear 
rate 7 varies from 10 -6 to 0.03, covering the linear and 
nonlinear response regimes of rhcology. For each state 
point under shear, 16 independent runs were performed 
to improve the statistics. 
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FIG. 1: The self part of (a) the intermediate scattering func- 
tion F 3 (q p ,t) and (b) the four-point susceptibility xVi't) at 
representative area fractions. 
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FIG. 2: ^-dependence of a-relaxation time r a , peak time of 
the four-point susceptibility r x and self-diffusion constant D 3 . 
The red solid lines are the results of the power-law fittings 
T a ~ {<t>c — <j>) ' with 7 = 2.7 for r a , 7 = 2.6 for r x and 
7 = 2.4 for D7 1 . Excellent fits are obtained to all data for 
a single value cj>c — 0.656 ± 0.001. The green dashed lines 
are the results of the Vogel-Fulcher-Tamman fittings r a ~ 
exp[D(j>/((j}o - <j>)] with D = 0.22 for r a , D = 0.23 for r x and 
D = 0.25 for D' 1 . Excellent fits are obtained to all data for 
a single value <j> = 0.671 ± 0.004. 



A. Slow dynamics and growing dynamic 
heterogeneity approaching the glass transition 



III. SLOW DYNAMICS AND GROWING 
LENGTH SCALES AT EQUILIBRIUM 

In this section, we first describe slow dynamics and 
growing dynamic heterogeneity approaching the glass 
transition point in a quiescent condition, and then dis- 
cuss growing dynamic and static length scales and struc- 
tural features of the slow dynamics. Relation between 
structure and dynamics in polydisperse colloidal liquids 
has also been discussed recently in Refs. [lj| and the 
emergence of long-lived clusters of particles with higher 
order can be connected to slow dynamics and dynamic 
heterogeneity Structural features of slow dynamics and 
dynamic heterogeneity are also confirmed by our results. 



We first consider the <j) dependence of the potential 
energy E p and the pressure p of the system. We find 
that both E p and p are smooth functions of <f> (data not 
shown), indicating that crystallization indeed does not 
occur. Further confirmations for the avoidance of crys- 
tallization come from the positional order and the spatial 
correlation of hexatic order, which will be shown later. 

The dynamic slowing down can be characterized by 
calculating the self part of the intermediate scattering 
function (ISF): 

1 N 

Fs(q P , t) = jj<Y, ex P^P ' ~ r,(0)]} >, (3) 

where < • • • > indicates the thermal average, i = y—T 
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and the wave number q p corresponds to the first peak of 
the static structure factor (which will be shown later). 
F s (q p , t) can reflect structural relaxation process of a liq- 
uid and representative results for several area fractions 
are shown in Fig. 1(a). We observe that the relaxation 
slows down and becomes more stretched as <j) increases. 
In particular, a two-step decay emerges at high enough 
area fractions, reflecting the rattling motion of parti- 
cles trapping within cages at short times (/3-rclaxation) 
and the motion of particles escaping from the cages at 
long times (a-rclaxation) . We define a- relaxation time 
as F s (q p ,t = r a ) = 0.2 (we have checked that choos- 
ing other values will not alter the qualitative results), 
which is plotted as a function of 4> in Fig. 2. Clearly, r a 
drastically increase as <p increases. For the density range 
investigated, r Q can be well fitted by the mode-coupling- 
theory (MCT) power-law [l?}: t q ~ (4> c — 0) -7 : where <f> c 
is the MCT glass transition point, or the Vogcl-Fulcher- 
Tamman (VFT) law: r Q ~ cxp [£></>/ (0o — <f>)], where D 
is the fragility index and <f>o is the ideal glass-transition 
point. 

For the glass-forming liquids, the dynamic slowing 
down accompanies growing dynamic heterogeneity. We 
quantify the dynamic heterogeneity by the self part of 
the four-point density correlations Q, which dominates 
results of the total four-point density correlations [H, [l7| ■ 
A time-dependent self-overlap order parameter Q s (t) is 
denned as 

1 N 

Q«(*) = jV<£™(M*)-*i(°)l)>> < 4 > 

with w = 1(0) for \rj(t) - r.,(0)| < (>)0.4 < a > . The 
threshold 0.4 was determined by considering the peak 
height of the four-point susceptibility Xa (*) = j^[< 
Q s (t) 2 > - < Q s (t) > 2 ] becoming maximal for = 0.64. 
X% s (t) measures the extent to which the dynamics at any 
two points in space are correlated within a time inter- 
val. Typical results for xT(t) are shown in Fig. 1(b). 
As can be seen, the peak height of xfW increases and 
the peak time r x shifts to larger times as <f> increases, 
which indicates that the dynamics not only slows down 
but also becomes progressively more heterogeneous on 
the approach to the glass transition. t x can also be well 
fitted by MCT power-law or VFT law, as shown in Fig. 
2. We also present <fi dependence of the self-diffusion 
constant D s in Fig. 2, as evaluated by using Einstein 
relation from the long-time data of the mean squared 
displacements. Independent fits give these three types 
of data with single parameters </> c = 0.656 ± 0.001 and 
</>o = 0.671 ± 0.004. 

As mentioned in Sec. I, the development of dynamic 
heterogeneity suggests the existence of a growing dy- 
namic length scale, which can be investigated by the 
four-point, time-dependent structure factor for the self- 
overlapping particles, which is defined as 

Sr(q,t) = ^<p(q,t)p(-q,t)>, (5) 
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FIG. 3: The structure factor of the self-overlapping particles 
Sl s {q,t) at low-q region for several area fractions. The solid 
lines are the results of the OZ fittings (see text). 
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FIG. 4: (a) The pair correlation function g(r) and (b) the 
static structure factor S(q) at low-<j region at varying cf). The 
results have been shifted for clarity. 

where p{q,t) = T,f =1 w(\rj(t) - r,-(0)|) exp[iq • ^(0)]. 
Here, the time t is often taken as r x , where dynamic 
heterogeneity becomes most pronounced Represen- 
tative results are given in Fig. 3. S'l s (q 1 T x ) at low-q 
region can be fitted by the Ornstein-Zernike (OZ) func- 
tion Sq/[1 + (g^) 2 ] and then we can obtain the dynamic 
correlation length £4 in this way. We will show the results 
of £4 later and discuss its link to the slow dynamics. 



B. Growing length scales and structural features of 
slow dynamics 

In Fig. 4, we present results of the pair correla- 
tion function g(r) and the static structure factor S(q) 
at varying <f>, as calculated by g(r) = 2 ^ rAr L N{N _ 1} < 
ZjjL k 5(r - \r jk \) > and S(q) = ^ < p(q)p(-q) > respec- 
tively, where p(q) = Y]j—i exp(iq ■ r^). On one hand, as 
found in other glass-forming liquids, g(r) and S(q) only 
display slight change on approaching the glass transition 
(notice that r Q increases by ~ 2 orders of magnitude in 
the studied density range, as shown in Fig. 2). On the 
other hand, although the long-range positional order is 



FIG. 5: (a) ^-dependence of *?6- The line is a linear fit to 
the data, (b) (^-dependence of ge(r)/g(r). The solid lines 
are the results of the OZ fittings of the envelopes. The green 
dashed line is included to stress the fact that crystallization is 
avoided since ge(r)/g(r) decays much faster than r~ 1,/4 . Note 
that ge(r)/g(r) decays as r~ 1//4 at the boundary between the 
liquid phase and the hexatic phase [5^ ]. 
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FIG. 6: (^-dependence of £4 and ^6- The red solid lines are 
the results of the MCT fittings ~ (0 C - with 4> c = 0.656 
and 7 = 0.83 for f 4 , </>c = 0.658 and 7 = 0.39 for f 6 . 
The green dashed lines are the results of the VFT fittings 
~ exp[_D0/(0 o - 4>)] with (j> = 0.670 and D = 0.074 for f 4 , 
0o = 0.678 and D = 0.048 for f 6 - The fitting results of £4 and 
£e for <f>c and 4>o agree well with that of r a , r x and D 3 within 
statistical accuracy. The dotted lines are the fittings of the 
functional form £ = £o[0/(0o — (/>)]■ 



obviously prevented from Fig. 4, the splitting second 
peaks in both g(r) and S{q) become more apparent as </> 
increases, suggesting the development of the locally pre- 
ferred order with increasing density. 

To characterize the local structure, we used a sixfold 
bond-orientation order parameter for each particle, which 
is denned as ipQ = ^- J2m=i ex P(*6^„). Here, rij is the 
number of the nearest neighbors for particle j which de- 
termined by the Voronoi construction, and 9 J m is the an- 
gle between (r m —Tj) and the x axis (particle m is a neigh- 
bor of particle j). By construction, 'J'g = |^g| is equal 
to 1 for a perfect hexagonal arrangement of six neighbors 
around particle j and for a random arrangement. Then, 
the time averaged order parameter over r Q for particle j is 
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FIG. 7: Log-log plot of £4 and ^6 versus r a . The lines are the 
results of the power-law (£ ~ r£) fittings with /j, — 0.28 for £4 
and fi = 0.11 for ^6- 



denoted by *g = £ J*° +T " ¥ 6 dt. We also calculated the 

order parameter of the system ^> 6 = j? < Ylj=i 
which is shown in Fig. 5(a) as a function of <f>. Obvi- 
ously, ^6 monotonically increases in the studied density 
range as <f> increases, indicating again the development 
of the locally preferred order on increasing density. Here 
we note that hexatic ordering in hard disc systems is a 
direct consequence of dense packing and manifestation of 
low configurational entropy, as emphasized in Refs. [l9| . 
The spatial correlation of ipg was then calculated as 

9s(r) = 2 ^Ar^(jV-i) < >■ The 

results of g§{r)/g{r) for four different area fractions are 
shown in Fig. 5(b). The envelops of g6{r)/g(r) can be fit- 
ted by the OZ function r -1 / 2 exp(— r/^), and the static 
correlation length can be obtained in this way. We 
show results of together with £4 (as estimated from 
Sl s (q,t)) in Fig. 6. Clearly, they both increase as cf> 
increases (or as the dynamics slows down). However, 
our data cannot be well described by the functional form 
£ = £o[</>/ (0o — 4*)} (see dotted lines in Fig. 6), which was 
proposed in Refs. [Ijj]. Instead, they can be well fitted 
by MCT or VFT functions as shown in Fig. 6, in good 
agreement with the results of Ref. [H, 0|. In fact, we 
find that the correlation length scales £ have power-law 
dependence on r a (i.e., £ ~ t£), as evidenced from Fig. 
7, where log-log plot of £4 and versus r a is presented. 
We will show later that this power-law relation £ ~ t£ 
holds even under shear. Thus, there are indeed structural 
features of the dynamic slowing down. 

To further understand the relation between the slow 
dynamics and the structure, we also calculated the relax- 
ation of particles with different structural order, as shown 
in Fig. 8(a). It is observed that the particles with high 
structural order dominate the relaxation process. Inter- 
estingly, the particles with 1 J? 6 G (0.2, 0.4) relax more 
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FIG. 8: (Left) Main: The self part of the intermediate scat- 
tering function F a (q p ,t) for particles with different structural 

order at (j> = 0.65. Here, is averaged over the time win- 
dow. Inset: the distribution of neighbors for particles at vary- 
ing \&g. (Right) Corresponding particle trajectory during an 
interval of T a . Different colors correspond to different struc- 
tural order averaged over an interval of r a . Clearly, particles 
with the highest structural order (magenta particles) are less 
mobile than the others and dominate the slow tail of F s (q p ,t). 



slowly than those with \l/g £ (0.4, 0.8), which is also con- 
firmed by the particle trajectory within a time interval 
of r a in Fig. 8(b), i.e., the cyan particles are less mobile 
than the green and yellow particles, suggesting the ex- 
istence of an unexpected relationship between hexagonal 
structure and dynamic heterogeneity. This indicates that 
particles with specific structures may responsible for the 
slow relaxation process, which has been pointed out re- 
cently in a study of binary soft disks [5j| . This feature is 
also evidenced by the inset of Fig. 8(a) where it reveals 
that the cyan particles mainly consist of particles with 
five and seven neighbors. 

Therefore, the structural features do exist in our model 
and the slow dynamics accompanies the development of 
heterogeneity and growing correlation length scales in the 
absence of shear. We focus on the correlation between the 
slow dynamics and the structure under shear flow in the 
next section. 



IV. CORRELATION BETWEEN 
NONEQUILIBRIUM DYNAMICS AND 
STRUCTURE UNDER SHEAR 

When a glass-forming liquid is sheared, marked shear- 
thinning behavior is expected to appear at sufficiently 
large shear rates. It is also found in the previous work 
that the shear viscosity 77 and T a change with shear rate 
7 with r) oc r a ~ 7^ (y > 0). However, it is not clear 
how correlation lengths vary with 7 and whether there 
exists a connection between dynamics and structure in 
shear flow. In this section, we first investigate the Theo- 
logical behavior and nonequilibrium dynamics of the sys- 
tem under steady shear and then discuss their structural 
features. 

We first present the time evolution of the sample- 
averaged shear stress < a xy > after startup shear in Fig. 
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FIG. 9: Time evolution of (a) the shear stress < a xy > and (b) 
the bond-orientation order parameter after startup shear 
for several shear rates at cj> = 0.65. The results are averaged 
over 200 independent runs. 
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FIG. 10: (a) 7-dependence of 77 at varying <f>. The arrow 
indicates that the shear-thinning behavior starts at a shear 
rate much slower than l/r Q at equilibrium (7 = 0). (b) The 
same data presented in the form of 77 versus (f> at varying 7. 
The solid lines indicate that the weak ^-dependence of 77 at 
high shear rates crosses over into the strong ^-dependence of 
77 (olive line: 77 ~ exp[0. 20/(0. 667 — </))]) a t l° w shear rates. 



9(a). Here, a xy is defined as 



>xy 



1 N 

7<£ 



PjxPjy 



N 



3=1 J>k 



dr 



(6) 



It is seen that the system reaches the steady state after 
a transient period (typically, it needs to be significantly 
larger than 7 _1 ). Before entering the steady state, a 
peak emerges at some yielding shear stress. This phe- 
nomenon is usually called overshoot behavior (or yield- 
ing behavior). Although it was commonly found in many 
systems [43|, [6(1 HJ , its detailed mechanism remains elu- 
sive. Here we show in Fig. 9(b) that overshoot behav- 
ior also accompanies structural change in our system. 
The viscosity of the system was then calculated from the 
steady state shear stress by using the constitutive equa- 
tion ?] = <T X y/j. The shear rate dependence of ri is given 
in Fig. 10(a). The same data are also represented in 
the alternative form of 77 versus <p in Fig. 10(b). Clearly, 
the rheological behavior becomes Newtonian at low shear 
rates where 77 approached its equilibrium value (7 = 0), 
and non-Newtonian at high shear rates where marked 
shear thinning (77 ~ behavior displays. It should be 
noted that the shear-thinning behavior starts at a shear 
rate much smaller than the inverse of a-relaxation time 
scale at equilibrium, which is consistent with previous re- 
sults [ll|, [43| and indicates that there may be a structural 
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FIG. 11: Shear rate dependence of (a) F s (q p ,t) and (b) xi"(t) 
at <t> = 0.65. 
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FIG. 12: (a) 7-dependence of r a at varying 0. (b) The same 
data presented in the form of r a versus (j) at varying 7. The 
solid lines indicate that the weak ^-dependence of r a at high 
shear rates crosses over into the strong (^-dependence of r a 
(olive line: r a ~ exp[0.24</>/(0.669 — (f>)]) at low shear rates. 
The larger symbols indicate the crossover area fractions for 
several shear rates which distinguish equilibrium dynamics 
from shear-controlled dynamics (see text). 
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FIG. 13: (a) 7-dependence of t x at varying <j). (b) The same 
data presented in the form of r x versus <f> at varying 7. The 
solid lines indicate that the weak (^-dependence of r x at high 
shear rates crosses over into the strong ^-dependence of r x 
(olive line: r x ~ exp[0. 250/(0. 671 — <j})]) a t l° w shear rates. 



relaxation process much slower than r a characterizing the 
decay of the two-body correlation. Meanwhile, we find 
that rj exhibits weak 0-depcndcnce at high shear rates 
and strong ^-dependence at low shear rates. 

In order to investigate the noncquilibrium dynamics 
under steady shear, Eqs. (3) and (4) are extended to 
the system under steady shear, which have the following 
form: 



F s {q,t) = — < ^exp{i[(q p -7iq pa .y)-rj(i)-q p -rj(0)]} >, 

(7) 



FIG. 14: (a) 7-dependence of £4 at varying (/>. (b) The same 
data presented in the form of £4 versus <f> at varying 7. The 
solid lines indicate that the weak (^-dependence of £4 at high 
shear rates crosses over into the strong (^-dependence of £4 
(olive line: £4 ~ exp[0.11<?!>/(0.676 — (j>)]) a t l° w shear rates. 
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FIG. 15: (a) Shear rate dependence of (a) g(r) and (b) g§{r) 
at d> = 0.65. 



and 



1 N 

s(t) = ^<Y d w(\r j (t)-v j (0) 



7 %(0)x|)>. (8) 



The results of F s (q p ,t) and xl s (t) for several shear rates 
at 4> = 0-65 are given in Fig. 11. Their overall shapes 
are very similar to those in quiescent conditions and the 
slow decay of FJq p , t) satisfies the time-shear superposi- 
tion property |38[ ■ We present the results of r a and r x in 
Figs. 12 and 13. Clearly, they share strong similarities 
with the viscosity: (1) t q a r x ~ 7 _!y also holds in the 
marked shear-thinning regime, although the exponent v 
is different for 77, r a and r x ; (2) 77, r Q and t x all exhibit 
weak (^-dependence at high 7 and strong <p dependence at 
low 7; (3) for sufficiently low 7, there exists a crossover 
area fraction 0(7) distinguishing equilibrium dynamics 
from shear-controlled dynamics (see larger symbols in 
Fig. 12(b)). Specifically, 77, r a and r x are nearly inde- 
pendent on shear rate below ^(7) and strongly affected 
by both <fi and 7 above 4>(j) ■ This phenomenon is in good 
agreement with findings of Ref . [62| . 

We can also obtain a dynamic correlation length £4 
from the four-point structure factor (Eq. (5)) [(33|, which 
is shown in Fig. 14. £4 also decreases in the marked shear 
thinning regime as 7 increases. However, contrary to the 
previous study (44j-|48|, we find that power laws fail to 
describe its dependence on 7 even in the marked shear 
thinning regime. 

Turing to the static properties, we present the shear 
rate dependence of g(r) and ge (r) in Fig. 15. We find 
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FIG. 16: (a) 7-dependence of £e at varying (f>. (b) The same 
data presented in the form of ^6 versus <f> at varying 7. The 
solid lines indicate that the weak ^-dependence of ^6 at high 
shear rates crosses into the strong (^-dependence of ^ (olive 
line: ^6 ~ exp[0. 050/(0. 678 — (j>)]) at low shear rates. 



plied forces overwhelm random thermal forces and the 
system exhibit behavior like a normal liquid (e.g., the 
relaxation process becomes fast and only weak dynamic 
heterogeneity displays under very strong shear, see Fig. 
11). However, the relation £ ~ t£ between £4, £e and r a 
retains for 7 <~ 10~ 3 , where both thermal fluctuations 
and applied forces are important in determining prop- 
erties of the liquids and the dynamics slows down and 
accompanies strong dynamic heterogeneity. Therefore, 
structure features for the slow dynamics can indeed exist 
even under shear flow. 



CONCLUSIONS 
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FIG. 17: Log-log plot of £4 and ^ versus r a at varying 7. 
The lines are the same as those in Fig. 7. 

that only weak changes display in the radial distribution 
function (Fig. 15(a)), while pronounced changes are ob- 
served in the bond-orientation correlation function (Fig. 
15(b)). ge{r)/g(r) under shear can also be well fitted 
by the OZ function and then we can calculate the static 
correlation length under shear, which is shown in Fig. 
16. Clearly, power laws also fail to describe 7 dependence 
of £ 6 . We also confirm that £4 and under shear have 
similar dependence on 7 although we didn't find an em- 
pirical or theoretical functional form for them. However, 
their connection to slow dynamics under steady shear is 
clear, as evidenced in Fig. 17, where we present log- 
log plot of £4 and versus t q at varying 7. As can be 
seen, the power-law relation £ ~ t£ only breaks down 
for very strong shear (7 >~ 10~ 2 ), where externally ap- 



In conclusion, we have investigated structural features 
of slow dynamics and dynamic heterogeneity in a 2D 
polydisperse glass-forming liquid in both quiescent and 
sheared conditions. We have confirmed the previous find- 
ings in the absence of shear that slow dynamics accom- 
panies the development of dynamic heterogeneity The 
correlation length scales £ grow with r Q as a power-law 
£ ~ t£. In the presence of shear, both r) and r a have 
power-law dependence on 7 in the marked shear thin- 
ning regime and exhibit weak ^-dependence at high shear 
rates and strong (^-dependence at low shear rates. The 
dynamic and static correlation lengths also decrease with 
increasing 7, but their dependence on 7 cannot be de- 
scribed by power laws in the same 7 regime. We fur- 
ther find that the relation £ ~ t£ between length scales 
and dynamics holds for not too strong shear where ther- 
mal fluctuations and external forces are both important 
in determining the properties of dense liquids. There- 
fore, our results demonstrate a link between slow dy- 
namics and structure in glass-forming liquids even under 
noncquilibrium conditions and can help to better under- 
stand the non-Newtonian behavior of supercooled liquids 
and glasses. 
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